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Numerical Stochastic Perturbation Theory. Convergence and features of 
the stochastic process. Computations at fixed (Landau) Gauge. 
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Concerning Numerical Stochastic Perturbation Theory, we discuss the convergence ol the stochastic process 
(idea ol the prool, leatures ol the limit distribution, rate ol convergence to equilibrium). Then we also discuss 
the expected fluctuations in the observables and give some idea to reduce them. In the end we show that also 
computation ol quantities at fixed (Landau) Gauge is now possible. 



1. Introduction 

Numerical Stochastic Perturbation Theory 
(NSPT) was introduced in ||l| and had successful 
applications [|| ||^ . Here we want to describe the 
main features of the underlying stochastic pro- 
cess. In the last part we show that a deeper un- 
derstanding of the process allows also computa- 
tions at fixed (Landau) gauge. 

NSPT was developed in the context of Stochas- 
tic Quantization Q. In this approach per- 
turbation theory is performed through a for- 
mal substitution of the expansion Urj{x^t) — > 
'^f^ g''U^''\x,t), in the Langevin equation: 



d_ 

di 



(V is a suitably defined group derivative). The 
Langevin equation in the algebra variables reads: 

^A^,(r;,x,0 = i?fF,^^(r;,x,t)+77^(x,0 

which gives a system of equations that can be 
solved numerically. 



2. Convergence and 
stochastic process. 



features of the 



We now consider P„[^("\ . . . , i] : the dis- 
tribution function of the first n perturbative com- 
ponents of the fields at Langevin time t. Pn cer- 
tainly exists for each n and t fixed. But we need 
to know: whether the limit distributions exist for 



t — > oo, to what kind of distribution they con- 
verge and at which rate. 

At this stage the answer to the first question is 
certainly not. This is due to two reasons which 
become clear if we look at the formal solution of 
the (perturbatively developed) Langevin equation 
in the Fourier transformed space (with zero initial 
conditions): 



4")(fc,t) = T^. / dse-'='(*-^)/(")(fc,5) 



(1) 



Colour indexes are left out, T^t^ and L^^ are the 
transverse and longitudinal abelian projectors, 

and/(")(fc,t)=5/f^"-'^(fc,t)V/^^^""'^(fc,t); 
fi°\k,t) = T]^{k,t). With /(■) we mean the 3 (4) 
gluons interaction terms, which only depend on 
the fields till the n—1 {n~2) perturbative order. 

This process must behave like a random walk 
in correspondence of both the gauge degrees of 
freedom and the mode fc = 0, since such degrees 
of freedom have no attractive force and conse- 
quently no damping factor. Thus they produce 
diverging fluctuations, even if their mean value 
shall be zero since the equivalence of Stochastic 
and Canonical quantization is true also when ex- 
panded in perturbation theory. 

The idea is to introduce an attractive force 
which keep the norms of the fields under con- 
trol, without affecting the observables ||]. To 
this end we interlace each step of Langevin dy- 



namic with a step of gauge transformation de- 
fined by: UJ^{x) = e™(=")C/^(a;)e-"'(^+'') (where 
w{x) = —^J2fj. ^/j^A') O'^s can prove that, by do- 
ing this, the system gains a force that drives it 
towards the Landau gauge (provided it is within 
the first Gribov horizon) [^. 

We will come back to this point later when con- 
sidering quantities in a fixed gauge. For the mo- 
ment suffice it to say that all gauge degree of free- 
dom now have an attractive force. 

The problem of zero modes appears also in 
usual lattice perturbation theory: the propagator 
is singular at fc = 0, and the common prescrip- 
tion is to exclude the zero-momentum degree of 
freedom of the fields. 

In our context the problem is a bit more subtle 
since we work in configuration space. Moreover 
if we set the fc = mode to zero in the gaussian 
noise r]'^{x,t), the free fields A^^^x,t) will have 
no = mode too, but this won't be true for 
the higher perturbative components. In fact the 
interaction introduces a zero mode contribution 
in the fields A^p^ {x, t) even if it was not present 
in the lower perturbative order. For instance the 
three gluons term gives: 

^M(0>*)l3si« = 2(27r)" J '^P^'i^(P + 'i) 

Ali~p,t)AU^q,t)v^l{0,p,q) 

{vlJa{k,p,q) = (5^t,(fc - + perm.). Wc must 
subtract these modes by hand at each step. 

In principle one should subtract the zero mode 
component after the updating of each single per- 
turbative component of each single link, and be- 
fore evaluating the next perturbative order. But 
this would be extremely expensive. It is conve- 
nient, instead, to subtract the zero mode after 
each sweep of the whole lattice. In this way we 
just introduce another error of order r, which is 
then extrapolated to zero. 

Once implemented the two corrections de- 
scribed above, it is possible to prove the conver- 
gence of the process. By that we mean that any 
correlation function of any perturbative compo- 
nent of the fields {(Yij has a finite 
limit for large Langevin time. 

We just give the main ideas and results. We 



first remark that all correlations of free fields 
converge at least like 0{e~'^ *) (if q is the lower 
momentum). Then we proceed by induction: It 
is convenient to re-write the solution (|l|) in dis- 
cretized Langevin time (t = Nt) distinguishing 
the memory of the past from the new contribu- 
tion of the random process: 

4")(fc,i) = e-''^Al"\k,t-T) + V^r^^{k,t) 
A^Hk,t) = e-'^'^A^\k,t~T)^Tfl^\k,t) 

The insertion of this formula into a correlation 
function reduces it into others of lower perturba- 
tive order. It is not difficult to evaluate the sum 
of the relevant part which survives in the limit 
T ^ Q aX t = Nt fixed, and then take the limit 
t —t oo. 

The result is the following: if every degree 
of freedom has an attractive force, as described 
above, then a limit distribution exists for each P„, 
and all their moments are finite (i.e. any corre- 
lation function of any perturbative component of 
fields is finite). Moreover convergence to equilib- 
rium is damped by a factor Ve~ where k is the 
lower momentum contributing and p is the global 
perturbative order of the correlation function. 

3. About fluctuations: status and ideas for 
improvement 

The fact that the limit distributions of these 
objects produce correlation functions which are 
all finite is not sufficient. We need to have an 
idea of how much such correlations can grow for 
a high number of fields or perturbative order. In 
fact some applications of NSPT need a knowl- 
edge of the perturbative coefficient with an ex- 
tremely high precision. Since this is a stochastic 
method the result is known with an error which 
is essentially given by the intrinsic fiuctuations of 
the correlation function one needs to calculate. 

To gain some insight in this problem it is use- 
ful to think of the process in terms of the un- 
derlying gaussian process •q{x^t). Any perturba- 
tive component O^^'-' of an observable O may be 
seen as a sum of correlation functions of 77's. In 
fact (0(rt) ^ Y.AV.^■■■V.M) and {(O^P)f) ~ 
Ect Y.ni'n^i ■ ■ ■ V^mV,ti ■ ■ ■ iI^m)- The cr's are some 



choices from all the possible ry's in the process and 
M is a number. Both depend on the theory, on 
the observable and on the perturbative order p. 
Since the fluctuations of 0*^^^ clearly depend on 
the number M and choices a, it would be impor- 
tant to be able to say something about them. 

It is quite easy to determine the number M for 
a particular Theory and observable. For instance 
in the Xcj)^ theory for O[0] — 0^ the relation be- 
tween the perturbative order p and the maximum 
number of correlated rfs M is M = 2p + 2, while 
for the plaquette in gauge theory the relation is 
simply M = p. 

This information is widely insufficient to de- 
termine the size of the fluctuations. There are 
other elements that strongly influence the size 
of the fluctuations, but we can be only quali- 
tative about them. Consider for instance Xcj)"^: 
the field interacts with itself in the same point. 

0(1) ^ (0(O))3^ ^(2) _ (0(O))20(1) _ (^(0))5^ 

have a strong contribution of correlations of the 
kind {r]^)- We expect strong fluctuations quite 
soon. The situation for a— model is similar. Con- 
sider instead gauge theories: Interaction is given 
by product of fields of different colours and direc- 
tions. Al,^'^^ - gP'"'A'l!:°'>At^°'> etc. This makes 
the previous phenomenon much less severe. 
Remark: If we want to study a fixed observable at 
a fixed perturbative order we do not really need 
the process 77 to be gaussian. We just need a 
fixed number of its correlations to be gaussian. 
Higer moments could be chosen to be lower than 
those of a gaussian distribution. This is achieved, 
for instance, if one exploits combinations of Dirac 
delta's: p{x) = '^jWj5{x — Xj) (see for the 
general solution). 

4. Computation of quantities at a fixed 
(Landau) Gauge. 

The convergence of each correlation function 
imply in particular that not only gauge invariant 
quantities are computable but also those which 
depend on the gauge. 

Actually the gauge fixing procedure which is 
realized here is not that of Faddeev-Popov 
(which is possible but more expensive to perform 
|9|), but that introduced by Zwanziger (plus cor- 



rections of the order of the lattice spacing and 
of t). In fact the interlaced gauge transforma- 
tions described above are equivalent to add to the 
Langevin equation a force 

XDfd^Al = -X^-^^^ + Xgr'^'^A^d^Al 

(where S'g_f[A] = | / dxld^A^Y). Corrections of 
the order of the lattice spacing are present, since 
the formula above is valid only in the continuum. 
A correction of order r is expected to come from 
the procedure of interlacing a Langevin step with 
a gauge transformation. 

Although this kind of gauge fixing is not that 
of Faddeev-Popov, The Landau choice of gauge 
can be reproduced. There are at least two ways 
of doing this. The first one is natural but maybe 
not efficient: one can perform the calculation at 
different value of the ratio — = A/r and then ex- 
trapolate for large ^ . The second method consists 
in performing many gauge transformations on a 
thermalized configuration. This should drive the 
system towards a stationary point, where - in fact 
- the Landau gauge condition is satisfied. 
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